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We consider quantum transport of spinless fermions in a ID lattice embedding an interacting region (two sites 
with inter-site repulsion U and inter-site hopping td, coupled to leads by hopping terms t c ). Using the numerical 
renormalization group for the particle-hole symmetric case, we study the quantum conductance g as a function 
of the inter-site hopping td- The interacting region, which is perfectly reflecting when td — > or td — > oo, 
becomes perfectly transmitting if td takes an intermediate value r(U, t c ) which defines the characteristic energy 
of this interacting model. When td < t c \/U, g is given by a universal function of the dimensionless ratio 
X — td/r. This universality characterizes the non-interacting regime where r = t 2 c , the perturbative regime 
(U < t^) where r can be obtained using Hartree-Fock theory, and the non-perturbative regime (U > i 2 ) where 
r is twice the characteristic temperature Tk of an orbital Kondo effect induced by the inversion symmetry. 
When td < T, the expression g(X) — A(X + X -1 ) -2 valid without interaction describes also the conductance 
in the presence of the interaction. To obtain those results, we map this spinless model onto an Anderson model 
with spins, where the quantum impurity is at the end point of a semi-infinite ID lead and where td plays the role 
of a magnetic field h. This allows us to describe g(td) using exact results obtained for the magnetization m(h) 
of the Anderson model at zero temperature. We expect this universal scaling to be valid also in models with 2D 
leads, and observable using 2D semi-conductor heterostructures and an interacting region made of two identical 
quantum dots with strong capacitive inter-dot coupling and connected via a tunable quantum point contact. 

PACS numbers: 71.10.-w,72.10.-d,73.23.-b 



Introduction 

In quantum transport theory, a nanosystem inside which the 
electrons do not interact has a zero temperature conductance 
which is given (in units of the conductance quantum e 2 /h for 
spin polarized electrons, 2e 2 /h with spin degeneracy) by 



g= \t ns (E F )\< 



(1) 



where \t ns (Ep)\ 2 is the probability for an electron at the 
Fermi energy Ep to be transmitted through the nanosystem. 
This Landauer-Biittiker formula can be extended 1 ' 2 to an inter- 
acting nanosystem, if it behaves as a non-interacting nanosys- 
tem with renormalized parameters as the temperature T — > 0. 
However, this effective non-interacting nanosystem does not 
describe only the interacting region, but depends also on the 
presence of scatterers which can be outside, in the attached 
leads. This non-local aspect of the effective transmission 
\t n s (Ep) | 2 is characteristic of nanosystems inside which elec- 
trons interact and has been studied in ID models 3,4 using the 
density matrix renormalization group (DMRG), and in ID 5 7 
and 2D models 8 using the Hartree-Fock (HF) approximation. 

In this work, we study how the effective transmission of 
a nanosystem with perfect leads is renormalized by local in- 
teractions acting inside the nanosystem, using the numerical 
renormalization group (NRG) algorithm 9-13 and an inversion- 
symmetric interacting model (ISIM). This model describes the 
scattering of spin-polarized electrons (spinless fermions) by 
an interacting region (two sites characterized by an inter-site 



hopping term td, coupling terms t c and an inter-site repulsion 
U). Our study is restricted to the symmetric case (i.e. the case 
where ISIM is invariant under particle-hole symmetry). 

Firstly, we prove that ISIM, which is perfectly reflecting 
when the inter-site hopping term td — > or td — » oo, ex- 
hibits a peak of perfect transmission for an intermediate value 
T(U,t c ) of td- This scale t(L/, t c ) defines very precisely the 
fundamental energy scale of ISIM. HF theory gives correctly 
this peak of perfect transmission when U < i 2 , but does not 
give it when U exceeds £ 2 , showing the existence of a non- 
perturbative regime where the use of the NRG algorithm is 
required. In this non-perturbative regime, 



r{U,t c )=2T K , 



(2) 



where Tk is the characteristic temperature of an orbital Kondo 
effect induced by the inversion symmetry. 

Secondly, we show that the zero-temperature conductance 
g is given by a universal function g(X) of the dimensionless 
coordinate X — td/r. This function g(X) is independent of 
the choice of t c and U as far as td < t c VU. When td < r, 



l{X) = A{X + X- 1 )- 



(3) 



When t < td < t c yU, the conductance g{X) can be de- 
scribed by another function 



g(x) 



2.02 



0.74 



]n(2.8X) 



(4) 
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which is related to an exact result obtained by Tsvelick and 
Wiegmann 14 ' 15 with Bethe-Ansatz for the magnetization of 
the Anderson model at zero temperature. When t& exceeds 
t c \fU, the interaction becomes irrelevant, r is given by its 
non-interacting value t% X = td/t 2 c and g(X) = A(X + 
X^ 1 )^ 2 again. The conductance g is one example of the 
physical properties of ISIM which are given by universal func- 
tions of td/r at zero temperature. The low-energy effective 
one-body excitations provide another example for which we 
show the corresponding universal curves. 

In order to obtain these universal functions, it is useful to 
notice that the inversion symmetry of ISIM gives rise to a 
pseudo-spin, allowing to exactly map this ID spinless model 
onto an Anderson model with spins where the inter-site hop- 
ping td plays the role of a magnetic field, the quantum im- 
purity being at the end point of a single semi-infinite chain. 
Therefore, the behavior of ISIM as a function of td is related 
to the behavior of the Anderson model as a function of an ap- 
plied magnetic field h. 

The paper is organized as follows: SectionQ]introduces uni- 
versal aspects which characterize the Anderson model and are 
relevant for quantum dots where spin or orbital Kondo effects 
occur. The studied spinless model with inversion symmetry 
(ISIM) is defined in section [TT] and mapped onto an Anderson 
model with magnetic field in section [III] A second transfor- 
mation is performed in section llVl based on the usual logarith- 
mic discretization of the energy band of the leads, to get the 
final model used for the NRG study. The study is restricted 
to the case with particle-hole symmetry in section [V] The 
low energy excitations are first considered as a function of the 
temperature in section [VI] When U > ttT, T — t 2 being the 
level width of the scattering region, the 3 fixed points [free 
orbital (FO), local moment (LM) and strong coupling (SC)] 
characterizing the Anderson model without field are recov- 
ered when td is small enough. As td increases, the LM fixed 
point characterizing the model above the Kondo temperature 
Tjc disappears. In section \VU\ we study the low energy ex- 
citations of ISIM as a function of td in the limit T — >• 0. We 
find that they can always be described by a set of effective 
one-body excitations, showing that a continuous line of free 
fermion fixed points goes from the SC limit of the Anderson 
model (td = 0) towards a new simple limit: the polarized or- 
bital (PO) fixed point (t d oo). Between the SC and PO 
fixed points, we show in section fVIHI that there is always an 
intermediate value r(U, t c ) of td for which ISIM is perfectly 
transmitting. r(U,t c ) defines the fundamental energy scale 
of ISIM. In section HXl a simple method for calculating g(td) 
from the effective one body excitations characterizing ISIM 
when T — > is introduced. Using this method, we give in 
section |IX] the main result of this work, i.e., if one uses the 
dimensionless ratio X — td/r, the physical properties (con- 
ductance or effective one-body spectra) are universal and in- 
dependent of t c and U as far as td < t c \/U . This universal 
regime is divided in a first regime where the system is not 
very far from the SC fixed point (td < t, subsection IX At 
and where g(X) — 4(X + X^ 1 )^ 2 , followed by a second 
regime where g(X) is given by another universal function 
(r < td < t c -JTJ, subsection IX Bt . In the equivalent An- 



derson model, this second regime is characterized by the oc- 
currence of a magnetic moment. When td > t e vU, the in- 
teraction U becomes irrelevant and g(X) = A(X + X^ 1 )^ 2 
with X = td/t 2 (subsection lXCI ). In section lXll we show that 
HF theory gives the values of g obtained from the NRG spec- 
tra and the scale t(U, t c ), if U does not exceed t 2 . In contrast, 
HF theory fails to give perfect transmission if U > t 2 , show- 
ing the existence of a non-perturbative regime for ISIM. To 
obtain r in the non-perturbative regime, we first revisit in sec- 
tion |XlI] a method giving g from the difference of occupation 
numbers between the even and odd orbitals of the nanosystem. 
This method based on Friedel sum rule (FSR) contains an as- 
sumption. If ISIM is near the SC fixed point (td < r, non- 
perturbative regime), this assumption turns out to be justified. 
In that case, g(td) can be obtained from the impurity magne- 
tization m(h) of the equivalent Anderson model with a mag- 
netic field h at zero temperature. Exact results giving m(h) 
are reviewed in section IXIIII for the Anderson model. Using 
those results, we show in section IXIVl that t(U, t c ) = 2Tk in 
the non-perturbative regime, Tk being the characteristic tem- 
perature of the orbital Kondo effect exhibited by ISIM when 
U > t 2 . Moreover, a fit inspired from the exact behavior of 
m(h) in the local moment regime is used for describing the 
universal function g(X) when r < td < t c VU. Eventually, 
we summarize in section IXVl the universal aspects obtained 
using a simple ID model, and we conjecture that they can be 
extended to 2D models and observed in 2D semi-conductor 
heterostructures, where the nanosystem would consist of two 
identical quantum dots coupled by a quantum point contact. 



I. ANDERSON MODEL, KONDO PHYSICS, QUANTUM 
DOTS AND UNIVERSALITY 

The Anderson model describes a single site with Hubbard 
interaction U coupled to a 3D bath of conduction electrons. 
This is one of the quantum impurity models 13 which were in- 
troduced to study the resistance minimum observed in met- 
als with magnetic impurities. The Kondo problem refers to 
the failure of perturbative techniques to describe this mini- 
mum. The solution of these models by the NRG algorithm, 
a non-perturbative technique 9-13 introduced by Wilson, is at 
the origin of the discovery of universal behaviors which can 
emerge from many-body effects. Without magnetic field h 
and with particle-hole symmetry 10 , the Anderson model maps 
onto the Kondo Hamiltonian if U > 7iT, T oc t 2 , being the 
impurity-level width. In that case, there is a non-perturbative 
regime where the temperature dependence of physical observ- 
ables such as the impurity susceptibility is given by univer- 
sal functions of T/Tk, Tk being the Kondo temperature. If 
U < 7iT, the impurity susceptibility can be obtained by per- 
turbation theory. Universality characterizes not only the be- 
havior of the Anderson model as a function of the tempera- 
ture T, but also its behavior at T = as a function of an 
applied magnetic field h. Using the Bethe-Ansatz, Tsvelick 
and Wiegmann 1415 have obtained for the magnetization m(h) 
a universal function of the dimensionless variable Ii/Tk when 
T — >• 0. 
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The possibility to design artificial magnetic impurities in 
nanoscale conductors has opened 1617 a second era for quan- 
tum impurity models. Measuring the conductance g of quan- 
tum dots created by electrostatic gates, in a 2D electron 
g as ifr i i8 or m car b on nanotubes 19 , one obtains values which 
can be on universal curves as functions of T/Tk if there is a 
Kondo effect. Moreover, quantum dots open the possibility to 
study the Kondo effect as a function of the coupling between 
the impurity and the continuum of conduction electrons, and 
not only as functions of the temperature and of the magnetic 
field. As pointed out in Ref. 20, this gives the opportunity 
to do the spectroscopy of the Kondo problem. Notably, the 
weak to strong coupling crossover can be studied by vary- 
ing gate voltages, when metallic gates are used for creating 
quantum dots. Kondo physics was first related to the anti- 
ferromagnetic coupling between a magnetic impurity and the 
spin of the host's conduction electrons. This is why Kondo 
physics was first expected and seen 1617 in quantum dots with 
odd numbers of electrons, weakly coupled to leads. However, 
it was realized that a localized electronic state coupled to a 
continuum can give rise to a large class of different Kondo 
effects, including the original spin- 1/2 Kondo effect, various 
orbital Kondo effects and the SU(4) Kondo effect occurring if 
a spin Kondo effect co-exists with an orbital Kondo effect. 

In this framework, ISIM is a model which can be used for 
describing the quantum conductance of spin polarized elec- 
trons in an inversion-symmetric double-dot setup with strong 
capacitive inter-dot coupling, as a function of the inter-dot 
hopping td- For such a setup, td could be varied by electro- 
static gates if the two dots are coupled by a quantum point 
contact, and our study describes the effect of this coupling 
upon the orbital Kondo effect induced by the inversion sym- 
metry. Eventually, universal aspects of many-body phenom- 
ena characterize not only equilibrium quantum transport, but 
also non-equilibrium quantum transport which occurs in the 
presence of a large source-drain bias V s d- Measures of the 
conductance 18 and of the current noise 19 of Kondo dots have 
recently confirmed the expected universality if one measures 
T or V s d in units of Tk- We describe here another uni- 
versal aspect of linear quantum transport, i.e., the quantum 
conductance of a spin polarized inversion symmetric double- 
dot setup should be a universal function of the dimensionless 
inter-dot hopping td/Tx when T — > 0. 

Kondo physics is also at the origin of spinless models, as the 
interacting resonant level model 21 (IRLM) which describes a 
resonant level (Vdcftd) coupled to two baths of spinless elec- 
trons via tunneling junctions and an interaction U between 
the level and the baths. IRLM, which is now used for study- 
ing non-equilibrium quantum transport 21 ^, is related to the 
Kondo model, the charge states rid = 0,1 playing the role of 
spin states. Both ISIM and IRLM are inversion symmetric. 
However, the Zeeman field acting on the impurity is played 
by the hopping term td for ISIM, and by the site energy Vd for 
IRLM. Therefore, ISIM does not transmit the electrons with- 
out "field", while IRLM does. Though we study in this work a 
finite density of particles, let us mention that the two-particle 
scattering problem has been solved 23 for ISIM. 
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FIG. 1: (Color online) Studied setup (ISIM) where spin polarized 
electrons (spinless fermions) can be scattered by a nanosystem made 
of the 2 red sites (energy Vg, inter-site repulsion U and an internal 
hopping td)- The nanosystem is embedded by coupling terms t c into 
a ID lattice (hopping term th). 



II. INVERSION-SYMMETRIC INTERACTING MODEL 

The ISIM model is sketched in Fig.Q]and consists of a ID 
tight binding lattice (hopping term th) where a finite density 
of spin polarized electrons (spinless fermions) can be scat- 
tered by a central region made of 2 sites of potential Vg, with 
an internal hopping term td, and 2 coupling terms t c . The dif- 
ficulty comes from the presence of a repulsion of strength U 
which acts if the two sites of the central region are occupied. 

The ISIM Hamiltonian reads: 



H — H n s + H r + Hi 



(5) 



where the Hamiltonian of the central region (the interacting 
nanosystem) is given by 

H ns = -t d (clc t + cjc ) + Vg (n + ni) + U rc. 7ii . (6) 

c x and c x are spinless fermion operators at site x and n x = 
c x c x . The coupling Hamiltonian between the nanosystem and 
the leads reads 

H c = -t c (cl lCQ + c\c 2 + H.c.) , (7) 

while the leads are described by an Hamiltonian 



Hi = -t h (cUx+i + H - c -) ■ 



(8) 



x— — oo 



where means that x = —1,0, 1 are omitted from the 



summation. 



III. EQUIVALENT ANDERSON MODEL WITH 
MAGNETIC FIELD 



Because of inversion symmetry, one can map ISIM onto 
a single semi-infinite ID lattice where the fermions have a 
pseudo-spin and where the double-site nanosystem becomes 
a single site with Hubbard repulsion U at the end point of 
a semi-infinite lattice. This equivalent Anderson model is 
sketched in Fig. [2] To show this mapping, we define the 
fermion operators 



"\x = (cl x+1 +< 



2. 



(cL +1 -4)/V2, 



(9) 
(10) 



4 



Un^n l ; 



v G ±t d 



FIG. 2: (Color online) Equivalent Anderson model: Electrons with a 
pseudo spin (f= even, J,= odd) are free to move on a semi-infinite 
chain (hopping term th) with a quantum impurity (red site) attached 
(hopping term t c ) to its end point. The inter-site ISIM interaction 
becomes a Hubbard interaction Un^n^ between impurity orbitals of 
different pseudo-spins. The impurity potential Vq has now a Zeeman 
term ±t d . 



v G ±t d 



where the /c-dependent hybridization 

V(k) = -i c V2/Vsinfc (15) 
yields an impurity level width at Ep which is given by 



T = -^smk F , (16) 
th 



nk,a = d\ IJ d k(T ande fc = -2t h cosk. 

One can see that ISIM is identical to an Anderson model 
with a local magnetic field td which acts on the impurity only 
and gives rise to the Zeeman terms ±td in Eq. (fTTT i. There- 
fore, in the limit td — > 0, ISIM must exhibit an orbital Kondo 
effect if the equivalent Anderson model can be reduced to a 
Kondo model. The fact that the impurity is not coupled to a 
3D bath of conduction electrons, but only to a single semi- 
infinite ID bath changes only the proportionality factor of the 
hybridization function. We underline that the dimensionality 
of the considered baths of conduction electrons does not play 
a significant role in Kondo physics, such that the results of this 
study should hold if one attaches 2D bars or 3D strips instead 
of ID leads to the same nanosystem. 



FIG. 3: (Color online) Corresponding NRG chain: The quantum im- 
purity (red site) is now coupled via an hopping term t' c (Eq. dl9t ) 
to a ID lattice, where the sites are labelled by n and describe con- 
duction electron excitations of length scale A™' 2 ^ 1 centered on the 
impurity. The successive sites are now coupled via hopping terms 
(Eq. d!8t ) which fall off as A~ n//2 . Used discretization parameter 
A = 2. 



which create a spinless fermion in an even/odd (c/o) combina- 
tion of the orbitals located at the sites x and — x+ 1 of the orig- 
inal infinite lattice, (or a fermion with pseudo-spin a — e/o 
in the transformed semi-infinite lattice). a c / o x are the corre- 
sponding annihilation operators. Expressing H ns in terms of 
these new operators, one gets 



H„ 



(V G -t d )n c + (V G +t d )n + Un c n , (11) 



where n a = a) a 1 a a 1 and where the pseudo-spin "e" ("o") is 
parallel (anti-parallel) to the "Zeeman field" td- In terms of 
the operators 



4, 



^^TsinWz-l))^ (12) 



creating a spinless fermion of pseudo-spin a and momentum 
k in the transformed semi-infinite ID-lead, the lead and the 
coupling Hamiltonians can be written as 



k,<T 



and 



H e = ^2v(k)(4 tl d kia + H.c), 



(13) 



(14) 



IV. CORRESPONDING NRG CHAIN 

ISIM can be studied using the NRG procedur e) 10 ' 12 ' 13 devel- 
oped by Wilson for the Anderson model after minor changes. 
First, we assume V(k) ~ V(kp = 7r/2) and, using stan- 
dard NRG procedure, we divide the conduction band of the 
electron bath into logarithmic sub-bands characterized by an 
index n and an energy width 



d n = A-"(l-A- 1 ). 



(17) 



Throughout this paper, we use the discretization parameter 
A = 2. Within each sub-band, we introduce a complete set 
of orthonormal functions tp np (e), and expand the lead oper- 
ators in this basis. Dropping the terms with p ^ and us- 
ing a Gram-Schmidt procedure, the original ID leads give 
rise to another semi-infinite chain with nearest-neighbor hop- 
ping terms, each site n representing now a conduction elec- 
tron excitation at a length scale A n / 2 kp 1 centered at the im- 
purity. In this transformed ID model shown in Fig. [3] and 
hereafter called the NRG chain, the impurity and the N — 1 
first sites form a finite chain of length N, which is described 
by the Hamiltonian Hn, the successive sites n and n + 1 be- 
ing coupled by hopping terms x n which decay exponentially 
as n — » oo and are given by: 



c n = A-/ 2 



(1 + A~ 1 )(1-A-"- 1 ) 
2y/(l - A- 2Ar -!)(l - A- 2Ar=f y 



(18) 



The impurity is coupled to the first site of the NRG chain by 
an hopping term 



t',= 



tc 



(87r 3 )V4 



, /A(A + i) 



(19) 
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Since the length N is related 10 to the temperature T by the 
relation 



k B T 



1 + A- 



_ A -(N-l)/2 



(20) 



AT can be interpreted as a logarithmic temperature scale (N oc 
— log T), the large values of N corresponding to temperatures 
T small compared to the bandwidth th- 

The NRG chain coupled to the impurity is iteratively di- 
agonalized and rescaled, the spectrum being truncated to the 
N s first states at each iteration (We use N s = 1024 in this 
study). The behavior of ISIM as T decreases can be obtained 
from the spectrum of Hm as N increases, the bandwidth of 
Hn being suitably rescaled at each step. A fixed point of the 
renormalization group (RG) flow corresponds to an interval of 
successive iterations N of the same parity, where the rescaled 
many -body excitations Ej(N) do not vary. The fixed point 
is therefore characterized inside this interval by two spectra, 
one characterizing the even values of N, the other the odd val- 
ues. If it is a free fermion fixed point, Ej = J2 a £a > tne e « 
being effective one-body excitations, and the interacting sys- 
tem behaves as a non-interacting system with renormalized 
parameters td and t c near the fixed point. Moreover, if one 
has free fermions when T — > 0, the conductance g can be 
directly extracted from the NRG spectrum. 



V. RESTRICTION TO THE SYMMETRIC CASE 

Using the NRG procedure, ISIM can be studied as a func- 
tion of T for arbitrary values of its 6 bare parameters U, Ep, 
Vg, td, t c and th- Hereafter, we take Ep = and Vg = 
— U /2. This choice makes ISIM invariant under particle-hole 
symmetry, with a uniform density (n x ) =1/2. Very often, the 
infinite bandwidth limit (th —> oo) is assumed in the theory 
of quantum impurities. This corresponds to magnetic alloys 
where the bandwidth of the conduction electrons is large com- 
pared to the other energy scales of the model. In this study, we 
take th = 1, which defines the energy scale and allows us to 
consider also mesoscopic regimes where the scales U, td, t c 
or Vq can exceed th- In that case, Eq. ( TToT l gives for the levels 
±td of the isolated non-interacting nanosystem a width 



T = tt 



(21) 



when the nanosystem is coupled to leads. Our motivation to 
restrict the study to the symmetric case is not justified by phys- 
ical considerations, but mainly for the sake of simplicity, re- 
stricting the RG flow into a space of 3 effective parameters 
(U, t c ,td) only. Doing so, we proceed as Krishna-murthy, 
Wilkins and Wilson for the Anderson model, studying first 
the symmetric case 10 before considering later the asymmetric 
case 11 and its characteristic valence-fluctuation regime. 



VI. ROLE OF THE TEMPERATURE T 

When td = 0, ISIM is an Anderson model which has the 
RG flow sketched in Fig. [4] for the particle-hole symmetric 







FO 





SC 



FIG. 4: Line of free fermion fixed points [thick solid line in the plane 
(td,t c )], characterizing ISIM when T — > as td increases from 
td = (SC fixed point) towards td — > oo (PO fixed point). The 
FO, LM and SC free fermion fixed points and the RG trajectories 13 
followed by ISIM as T decreases for td — are indicated in the 
plane (U, t c ), for tiT > U (dashed) and tiT < U (solid). 
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FIG. 5: (Color online) Many body excitations Ei as a function of N 
(even values) for U = 0.005 and t c = 0.01. For t d = (Fig.EJl), 
one can see the 3 successive plateaus (FO, LM and SC fixed points) 
of the Anderson model. As td increases (Fig. [5J5 and Fig. [5};), the 
LM plateau shrinks and disappears when > t*d = t c yU. 



case. At low values of N (high values of T), ISIM is located 
in the vicinity of the unstable free-orbital (FO) fixed point. 
As N increases (T decreases), ISIM flows towards the stable 
strong-coupling (SC) fixed point. 

If the interaction is weak (U < irT), its effects can be de- 
scribed by perturbation theory, the flow goes directly from the 
FO fixed point towards the SC fixed point, and there appears 
no orbital Kondo effect for td — > 0. 

If the interaction is large (U > irT), the flow can visit 
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td = Strong coupling 



•-O 




Singlet state Free fermions 



td = oo Polarized orbital 




Free fermions 



FIG. 6: (Color online) NRG chain in the SC limit (upper figure) and 
in the PO limit (lower figure). When td — > 0, the impurity (red 
dot) and the first site of the NRG chain form a system in its singlet 
ground state decoupled from the other sites which carry free fermion 
excitations. When td — > oo, the even (odd) orbital of the impurity is 
occupied (empty) and the other sites carry free fermion excitations. 
Therefore, there is a permutation of the parity of the length of the 
free part as td increases. When the excitations of the free part are 
independent of this parity, td = t and g — 1. 
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FIG. 7: (Color online) One body excitations e a (td) [extracted from 
the E^N -> oo, td)] for U = 0.1 and t c = 0.1 (left scale). The 
solid (dashed) lines correspond to NRG chains of even (odd) length 
N. Conductance g(td) extracted from Ae(td) using Eq. J27b (thick 
red curve, right scale). For td = t, the e a are independent of the 
parity of N and g = 1. 



an intermediate unstable fixed point — the local-moment (LM) 
fixed point — before reaching the SC fixed point. In that case, 
ISIM is identical to a Kondo model characterized by a temper- 
ature Tk and by universal functions of the ratio T /Tk- For 
td = 0, ISIM is on the FO fixed point when t c \f~U < T, ex- 
hibits a local moment when Tk < T < t c \/\J and reaches the 
SC fixed point when T < Tk- While the Hartree-Fock theory 
qualitatively describes 24 the local moment at high tempera- 
tures Tk < T < t c \/U, it breaks down at low temperatures 
(T < Tk), where the effect of the interaction becomes non- 
perturbative and gives an orbital Kondo effect. 

In Fig. [5] the first many-body excitations Ej of ISIM are 
given for increasing even values of N for td = 0. Since 
U > irt 2 c , one gets 3 plateaus corresponding to the 3 expected 
fixed points. Inside the plateaus, the spectra are free fermion 
spectra which are described in Ref.[l0|. However, between the 
plateaus, there are no free fermion spectra and Ej ^ ^2 a e a . 
As td increases (Fig. |5), the LM plateau decreases and van- 
ishes when td reaches a value s» t c \U. 



VII. ROLE OF THE INTERNAL HOPPING AT T = 

Let us study how the many-body levels Ej given by the 
NRG algorithm depend on the internal hopping term td in the 
limit where N — > oo, i.e. in the limit where the temperature 
T ^ 0. 

When td = 0, one has the SC limit 10 of the Anderson model 
where the impurity is strongly coupled to its first neighbor in 
the NRG chain (the conduction-electron state at the impurity 
site). The impurity and this site form a system which can be 
reduced to its ground state (a singlet), the N — 2 other sites 
carrying free fermion excitations e a which are independent of 
that interacting system. This SC limit of the Anderson model 
without field is sketched in Fig. [6] (upper part). In the Kondo 



model, the site directly coupled to the impurity describes the 
cloud of conduction electrons which fully screens the mag- 
netic moment of the impurity. 

When td — > oo, the impurity occupation numbers (n e ) — > 1 
and (n Q ) — > 0, and the — 1 other sites of the NRG chain are 
independent of the impurity. This limit is sketched in Fig. [6] 
(lower part). We call this limit "Polarized Orbital (PO)," since 
it coincides with the FO fixed point of the Anderson model, 
except that the spin of the free orbital is not free, but fully 
polarized in our case. 

The Ej(td) correspond to many-body excitations of effec- 
tive non-interacting spectra when td and td — > oo. When 
td varies between those 2 limits, the NRG algorithm contin- 
ues to give many -body excitations Ej(td) compatible with the 
free fermion rule Ei(td) = J2 a e a(td), allowing us to ex- 
tract one-body excitations e Q (t t j) for intermediate values of 
td- We conclude that there is a continuum of effective non- 
interacting spectra which describe the Ej as td varies in the 
limit N — > oo, i.e. the T — !• limit of ISIM is given by 
a continuum line of free fermion fixed points. This line is 
sketched in Fig. [4] Having always free fermions as td varies 
means that the T = scattering properties of an interacting 
region embedded inside an infinite non-interacting lattice are 
those of an effective non-interacting system with renormal- 
ized parameters, in agreement with the DMRG study of the 
persistent current given in Ref. 2. We underline that those ef- 
fective non-interacting spectra describe the T = limit, while 
a description of the low-temperature dependence of the con- 
ductance requires effective Hamiltonians of Landau quasipar- 
ticles with residual quasiparticle interactions. Such a Fermi 
liquid theory has been proposed by Nozieres. In the case of 
the Anderson model, it has been developed in Ref. 25 without 
magnetic field (td = 0) and in Ref. 26 with magnetic field 
(td 0). However, the quasiparticle interaction comes into 
play at finite temperatures only and a residual interaction is 
not necessary for describing the Ei as a function of td in the 
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limit TV — > oo. 

Fig. [7] shows these first one-body excitations e a as a func- 
tion of td extracted from the E](td), calculated with U = 0.1 
and t c = 0.1. The pseudo-spin degeneracy being broken by 
the "magnetic field" td ^ 0, the first (second) one-body ex- 
citation ei (€2) carries respectively an even (odd) pseudo-spin 
if N is even. This is the inverse if N is odd, e\ (£2) carrying 
respectively an odd (even) pseudo-spin. 

VIII. PERFECT TRANSMISSION AND 
CHARACTERISTIC ENERGY SCALE 

Since the free part of the NRG chain has N — 2 sites for 
N -> 00 and t d -> (SC fixed point), while it has N — 1 
sites for t d — » 00 (PO fixed point), there is a permutation of 
the e a (td) as increases: the e a (td — > 0) for N even become 
the e a (td — > 00) for N odd and vice-versa. This permutation 
is shown in Fig. [7] Since for N — > 00 there is a permuta- 
tion between the even and odd spectra as td increases, there 
is a value of td for which the e a (td) are independent of the 
parity of N. This value defines very precisely the character- 
istic energy scale t(U, t c ) of ISIM. Because of particle-hole 
symmetry, the nanosystem (the impurity of the NRG chain) 
is always occupied by one electron. Binding one electron of 
the leads with this electron reduces the energy when td < r, 
while it increases the energy when td > T. For td = t, it is 
indifferent to bind or not an electron of the lead with the one 
of the nanosystem, making ISIM perfectly transparent. This 
gives the proof that, for all values of U and t c , there is al- 
ways a value t of td for which the interacting region becomes 
perfectly transmitting and 

g(t d = r(U,t c )) = l. (22) 

The argument is reminiscent to that giving the condition for 
having a perfectly transparent quantum dot in the Coulomb 
blockade regime: td in our case, the gate voltage in the other 
case, have to be adjusted to values for which it costs the same 
energy to put an extra electron outside or inside the dot. The 
quantum conductance g can be extracted from the NRG spec- 
tra. Using a method explained in the following section, we 
have calculated g(t d ). The result shown in Fig.|7]confirms that 
g = 1 precisely for the value t of td for which the N — » 00 
low energy excitations are independent of the parity of TV. 

IX. FREE-FERMION SPECTRA AND QUANTUM 
CONDUCTANCE 

As pointed out in previous works 27-29 , the quantum con- 
ductance g can be directly extracted from the NRG spectra. 
Let us consider an NRG chain of even length N — > 00. When 
td = 0, the one body spectrum is identical for the two pseudo- 
spins. A hopping td ^ breaks the pseudo-spin degeneracy 
and opens a gap Ae(td) (indicated in Fig. O between excita- 
tions of opposite pseudo-spins. For free fermions in the limit 
T — >• 0, the asymptotic electron states of ISIM are station- 
ary waves with even (odd) phase shifts S e>0 (k) induced by the 
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FIG. 8: (Color online) e a and g as a function of td/t^, without in- 
teraction (U = 0). The values of g extracted from the NRG spectra 
(red cross) coincide with the exact values (Eq. d28t - red line). 

scattering region 

* e (M oc cos(fc(p - 1/2) - 6 c (k)), (23) 
*o(M oc sin(A:(p - 1/2) - o (fc)) . (24) 

The shifts S 0yO (k) of the scattering phases and Se Ct0 (k) of the 
energy levels are proportional for each pseudo-spin. This can 
be shown by taking a finite size L for ISIM, quantizing the 
momenta fc c o (n) = irn/L + d ej0 (k)/L and using the disper- 
sion relation e(fc) = 2cos(fc). This yields 

Ae = (Se e — Se ) oc (<5fc c — 8k ) cx (S c — 5 Q ) . (25) 

In the limit L — > 00, the quantum conductance g(td) of ISIM 
can be expressed as a function of the scattering phase shifts: 

g{t d ) = sin 2 {S c (t d ) - 6 (t d )) . (26) 

From Eq. ( |25T ) and Eq. d26i i. one eventually obtains the relation 
which allows us to extract g from the NRG spectra: 

g(td) = sin 7T — - . (27) 
V Ae(i d ^oo)/ 

The proportionality factor between Ae(td) and S e (t d ) — 5 (td) 
has been determined from the condition that g —> (for 
<5c — <5o — > 7r) when td — > 00. Eq. j27l ) describes a quantum 
conductance which vanishes when t& — > and when td — > 00, 
and reaches the unitary limit for an intermediate value t(U, t c ) 
of t d . 

Let us check that the NRG algorithm and Eq. d27l ) give 
us the correct behavior for the conductance g in the non- 
interacting limit where it is straightforward to solve the scat- 
tering problem. One obtains 

The e a (td) given by the NRG algorithm for U = are shown 
in Fig. |8] with the corresponding values of g obtained from 
Eq. d27| i. One can see that the behavior of g given by the 
NRG algorithm and Eq. d27| i reproduces the correct behavior 
given by Eq. (f28b in the non-interacting limit. 
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FIG. 9: (Color online) g(td) extracted from the NRG spectra for 
t c = 1 and many values of U. Increasing U shifts the transmission 
peaks to smaller values of td- The curves correspond respectively to 
U = (red curve), 5, 10, 15, 20, 25, 30, 35. 



X. UNIVERSAL PROPERTIES OF ISIM 

In Fig. [9] the conductance g extracted from the NRG spec- 
tra using Eq. d27l ) is given as a function of td for a coupling 
term t c = 1 and many values of U. The larger is U, the 
smaller is the characteristic scale t(U). Fig. [9] seems to in- 
dicate that the left sides of the transmission peaks are simply 
translated to lower values of td as U increases. This is con- 
firmed in Fig. [TOh ) where g is given as a function of the di- 
mensionless scale td/r, r being obtained from the criterion 
g(td = r) = 1. The curves g(td/r) obtained for t c = 1 are 
shown in Fig. [TOk). while the curves obtained for t c = 0.1 
and t c — 0.01 have been added in Fig. [TUfb). These fig- 
ures show the main result of this study, i.e., When td is not 
too large, g is given by a universal function of td/T(U,t c ). 
This function is independent of the values taken for U or t c . 
Since g(td) has been directly extracted from the free fermion 
NRG spectra, the NRG spectra must be also given by univer- 
sal functions of td/r(U,t c ), independent of the values of U 
and t c . This is shown in Fig. fTOl c) (t c — 1) and in Fig. fTOl d) 
(t c = 1, 0.1, 0.01), where the first excitations e a obtained for 
NRG chains of large even length N are plotted as a function 
of td/r. 

The values of t(U, t c ) used in Fig.[l0]are given in Fig.fTTI 
If one decreases the temperature T in the Anderson model 
without magnetic field, it has been shown in Ref. 10 that the 
interaction effects remain perturbative when U < ttT while 
they become non-perturbative when U > ttT. Since the level 
width T — t 2 c for the nanosystem used in ISIM at half-filling, 
we give the dimensionless scale r/U as a function of U /t 2 . in 
Fig-Qj] One can see that t/U has a slow decay followed by 
a faster decay, with a crossover around an interaction thresh- 
old consistent with the interaction threshold irT characteriz- 
ing the perturbative-non-perturbative crossover in the Ander- 
son model. For ISIM, this suggests that the interaction effects 
upon g are perturbative when U < ttT, and non-perturbative 
when U > ttT. However, the universal behavior of g(td/r) 



is not restricted to the non-perturbative "Kondo" regime, but 
characterizes also the perturbative regime. 



A. Universality near the SC limit (td < t) 

The universal regime can be divided into two parts. The 
first one begins at td = in the vicinity of the SC fixed point 
and ends at td = t. In this SC regime, g behaves as without in- 
teraction, but with a renormalized level width (t 2 — > r{U,t c )). 
As shown in Fig.[l0] all the curves g(td/r) are on a single uni- 
versal curve when td < r, independent of the values of U and 
t c . Since one of those curves (the red one) corresponds to the 
non-interacting limit U = 0, the universal curve for td < t is 
given by Eq. (l28T > where td is measured in units of r, instead 
of# 



B. Universality around a LM limit (r < td < t c y/U) 

The universal regime persists when td exceeds r. While 
g(td/r) decays immediately after the transmission peak if 
[7 = 0, g(td/r) begins to follow a new part of the univer- 
sal curve if U 7^ 0. This new part is not given by Eq. j29l , 
and ceases when a faster decay occurs. The larger is U, the 
larger is the interval of values of td/r where g(td/r) follows 
the slow decay of the universal curve (see Fig. ITOb . To de- 
scribe this slow decay, one can use exact results which give 
the magnetization m(h) of the Anderson model at T = as 
a function of the magnetic field h. This will be done after a 
study of the relation between g{td) and m(h). 

Let us just note now that the singlet state of the SC limit 
could be broken either if the temperature T or the "Zeeman 
energy" td exceeds the Kondo temperature Tr. This makes 
likely that the effects of T and td would be somewhat similar. 
If this is the case, the intermediate values of td would be re- 
lated to the formation of a local moment in the equivalent An- 
derson model of ISIM, and an intermediate LM regime would 
take place between the SC regime for low values of td and the 
PO regime for large values of td- This classification is used 
in Ref. 14 for describing the effect of a magnetic field in the 
zero-temperature limit of the Anderson model. However, we 
can only refer to the SC, LM and PO regimes, and not to the 
SC, LM and PO fixed points. The RG flows of ISIM yielded 
by increasing td at T — and by increasing T at td = 
are very different. Increasing T in ISIM at td = yields the 3 
plateaus shown in Fig. [5] characteristic of 3 well-defined fixed 
points. There are no free fermions between the plateaus, since 
there are no e a (T) such that Ei(T) = J2 a e a{T) outside the 
plateaus. In contrast, Fig. [10] does not exhibit plateaus and 
the Eiitd) can be described by a continuum of effective non- 
interacting spectra as td varies at T — 0, and not only by the 
three spectra of the SC, LM and FO fixed points. 
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FIG. 10: (Color online) Figure [Tol a): g as a function of td/r for t c = 1 and many values of U. r has been determined from the criterion 
g(td = t) = 1. The larger is [/, the larger are the values of td/r where g decays. The curves correspond respectively to U — (red curves) 
and [7 = 5, 10, 15, 20, 25, 30, 35 (black curves). Figure [Mb): 

g as a function of td/r. To the data (black curves) calculated taking t c — 1 
and shown in Fig.QIjta), we have added the data calculated taking t c = 0.1 and U = (red curve), 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4 
(green dashed curves) and taking t c = 0.01 and U = (red curve), 0.0005, 0.001, 0.0015, 0.002, 0.0025, 0.003, 0.0035 (blue dotted curves). 
Figure [T0lc): First one body excitations e a for iV even as a function of td/r. Same values of t c and U as in Fig.llOfa). Figure [TOld): e a for 
N even as a function of td/r. Same values of t c and U as in Fig. llOf b). 



Interaction-independent conductance in the PO limit 

(t d > t c VU) 



haves as without interaction when td exceeds a threshold value 



As one increases td/r, g(td/r) eventually exhibits a fast 
decay which is not given by a universal function of td/r. As 
can be seen in Fig.[l0] this fast decay corresponds to the decay 
of g{td/t 2 c ) obtained without interaction, but shifted to val- 
ues of td/r which increase when U increases. According to 
Ref.[l4|, the Anderson model at T = is in a LM regime if the 
magnetic field h lies in the interval Tk < h < h* oc VLT. 
The upper threshold ^/UT appears in the Hartree-Fock study 
made by Anderson 24 of the transition from the non-magnetic 
to the LM regime of the Anderson model, as one decreases 
T without magnetic field. For ISIM, this suggests that the 
slow universal decay of g(td/ T ) corresponding to the LM 



regime persists as far as td < t j oc t c yU. Above f^, ISIM 
should enter in the PO regime. This is confirmed in Fig. [T2l 
where one can see that g becomes independent of U and be- 



XI. CHARACTERISTIC ENERGY SCALE IN THE 
PERTURBATIVE REGIME 

Without interaction, Eq. d28l l implies that g = 1 if td = 
t\. This yields for the characteristic energy scale t of ISIM 
a non-interacting value £ 2 . For weak values of U, there is 
a perturbative regime where g and r can be obtained using 
self-consistent Hartree-Fock theory. In the symmetric case, 
the Hartree corrections and the site potentials Vg cancel each 
others since Vg — —U/2. The value of the inter-site hopping 
td is modified because of exchange and takes 5 a value given 
by the self-consistent solution of the HF equation 



= td + U(clc 1 {v,t c )) 



(30) 



Using Eq. (128b with v instead of td gives the HF value of g. 
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FIG. 11: (Color online) Characteristic scale r(U,t c ) as a function 
of U and t c - The values of T(U,t c )/U obtained from the condition 
g(td/r) — 1 and used in Fig. [TO] are plotted as a function of x = 
U/tl (+)• The solid blue line y(x) = Q.72%^2/ (-kx) exp -(ttx/8) 
fits the data in the non-perturbative regime and corresponds to the 
relation r = 2T K with T k = 0.364^2^(7/^ exp -(irU/Stl). 
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FIG. 13: (Color online) Conductance g as a function of id for tc — 1 
and {/ = 0, 0.5, 1, 1.5, 2, 2.5, 3. gr extracted from the NRG spec- 
tra (solid red curves) coincides with the HF estimates (dashed green 
curves) when U < or for i d > t e \/[7 if U > t\. 
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FIG. 12: (Color online) Ratio g(U,t c ,t d )/g(U = 0,t c ,t d ) as a 
function of t d / (t c y/U). g{U) is extracted from the NRG spectrum 
and g(U = 0) is given by Eq. l |28t . One can see that g behaves as 
without interaction when t d > t* d « lQt c ^/U . Same values of t c and 
U as in Fig. [H 



For the Anderson model, it is well known that HF theory 
fails to describe the Kondo regime. This Kondo regime oc- 
curs when the interaction exceeds a threshold value 7rr, ei- 
ther for low temperatures T < Tk without magnetic field 
h, or for weak fields h < Tk at T = 0. In the frame of 
the HF approximation, a magnetic moment should be formed, 
while it actually vanishes because of strong correlations be- 
tween the conduction electrons and the impurity spin (Kondo 
effect). Therefore, one does not expect that HF theory should 
be valid in the orbital Kondo regime of ISIM for large values 
of U (U > ttT) and small values of t d if T = 0. 

For a large coupling t c — 1, the results shown in Fig. Qj] 
confirm this prediction: HF theory gives the correct value of 
g for all values of td as far as U remains smaller than irij.. 
For U > Tit 2 c , the HF curves and the NRG curves coincide 
only when td is large (td > r), but become very different 
at the left side of the transmission peak. For a small coupling 



t c = 0.1, one can see in Fig.[l4]a more dramatic breakdown of 
HF theory which fails to give the peak of perfect transmission 
(see the curves with U = 0.04 and more notably U = 0.05 of 
Fig. [Hi. 

The origin of this failure can be simply explained. In the 
perturbative regime where g is given by Eq. d28l l with v in- 
stead of td, g = 1 if v = t;?. This yields for the scale r a HF 
value 

T HF =t 2 c -A(t c )U, (31) 

where the function 

A(t c ) = (cl Cl (v = t 2 c ,t c )), (32) 

which is shown in Fig. [15] depends weakly on t c (A = 1/tt 
if t c = 1 while A 1/4 if t c 0). When U reaches 
a value t 2 c /A, HF theory predicts that the interaction should 
renormalize td to a value thf = for having g = 1! This 
is absurd and confirms that HF theory breaks down for ISIM 
above an interaction threshold which is essentially the same 
as for the Anderson model (U > irT). 

In the non-perturbative regime of the Anderson model, the 
physical quantities such as the magnetization rn must be uni- 
versal functions of T /Tk or K/Tk- This means that t should 
be related to the characteristic temperature Tk of the orbital 
Kondo effect yielded by the inversion symmetry of ISIM. The 
relation between r and Tk can be obtained using Friedel sum 
rule and analytical results for the magnetization m(h) of the 
Anderson model at zero temperature. 

XII. IMPURITY OCCUPATION NUMBERS AND 
SCATTERING PHASE SHIFTS 

For having g, another method consists in using NRG for 
calculating the average impurity occupation numbers (n / ) 
of the even and odd orbitals of the nanosystem. The difference 



11 




FIG. 14: (Color online) Conductance g as a function of td for tc — 
0.1 and U = 0,0.01,0.02,0.03,0.04 and U = 0.05. The NRG 
results (solid red curves) and the HF results (dashed green curves) 
coincide for U < t 2 c . For U = 0.05, the HF curve gives only a small 
peak where g w 0.25, and not g = 1. 
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FIG. 15: Function A(t c ) defined in Eq. i32i as a function of t c . 




FIG. 16: (Color online) Values of g obtained from the NRG spec- 
trum (black solid lines, Eq. d27| >) and approximate values g (red 
dashed lines, Eq. d33t). as a function of td- Increasing U moves 
the conductance peak towards the left side. Figure (a): t c — 0.1 and 
[7 = 0, 0.01, 0.02, . . . , 1. Figure (b): t c = 1 and U = 0, 1, 2, 5, 10. 
Figure (c): t c = 0.1 and U = 0,0.05,0.1,0.15,0.2,0.25. Figure 
(d): t c = 1 and U = 0, 5, 10, 15, ... , 55. 



between the scattering phase shifts can be given in terms of 
(n / o ), if one assumes an approximate form of the Friedel 
sum rule (FSR): 

5 e - S fa 7r ((n c ) - (n Q )) , (33) 

and one gets an approximate value g for g from this estimate 
of S e — S Q using Eq. d26b : 

.g-sin 2 (7r((n c )-(n )) . (34) 

This approximate FSR is often used (see for instance 
Ref . 1301) and allows us to obtain g from the zero temperature 
impurity magnetization m(h) of the Anderson model with 
magnetic field h. This is particularly interesting since m(h) is 
a physical quantity for which exact results have been obtained 
with the Bethe-Ansatz by Tsvelick and Wiegmann. Unfortu- 
nately, Eq ( 1331 is only an approximation, and not the true FSR, 
as pointed out by Simon and Affleck in a study 31 of persistent 
currents through a quantum dot at Kondo resonance. The gen- 
eralization of FSR by Langreth 13 ' 32 shows that the phase shifts 
are proportional to the number of electrons displaced by the 
impurity, "among which are included not only the d electrons, 
but also some of the conduction electrons". For ISIM, this 
means that the displaced electrons are not only those inside 
the interacting region, but displaced electrons in the neighbor- 
ing parts of the leads have to be included too for obtaining the 
phase shifts from the occupation numbers via FSR. At first 
sight, one can expect that Eq. d33l could be used only if the 
scattering region is weakly coupled to the attached leads. Our 
results show that this is less simple. 

The difference between the values of g obtained directly 
from the NRG spectra [Eq. d27l il and the approximated values 
g are given as a function of td for weak (i c = 0.1, Fig.fTBh and 
d) and large values (t c = 1, Fig. [16b and c) of the coupling, 
and for weak (Fig.[T6h and c) or large values (Fig. [16b and d) 
of U. Even for a weak coupling t c = 0.1, where one could 
expect a negligible displaced charge outside the nanosystem, 
g = g only when ISIM is near the SC fixed point (td < r). 
When td > t, j / g. For a larger coupling (t c — 1), one 
can notice also differences between g and g even when td < t 
in the perturbative regime (U < iri^.). One concludes that g 
can be obtained from g with a good accuracy only in the non- 
perturbative regime where ISIM exhibits an orbital Kondo ef- 
fect (U > t 2 c jA and td < t). Otherwise, the difference S c — 6 
is not given by the difference (n c ) — (n Q ) evaluated inside the 
nanosystem, but depends also on the occupation numbers out- 
side the nanosystem. 

The validity of the approximate FSR in the non-perturbative 
regime when td < t can be explained by the following argu- 
ment: the conduction electrons which are displaced to screen 
the impurity pseudo-spin (forming a singlet state with the im- 
purity) are only a negligible fraction ps Tk/Ef of the con- 
duction electrons. To neglect this fraction induces an error 
oc T/th which cannot be seen in the curves shown in Fig. [TBI 
when t is very small (see Fig.[TTI). When the system is not in 
the orbital Kondo regime (U < irt 2 c or td > t), the number of 
displaced electrons becomes much larger, and the difference 
between g and g can be seen in Fig. [TBI 
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XIII. ZERO TEMPERATURE MAGNETIZATION OF THE 
SYMMETRIC ANDERSON MODEL 



XIV. CHARACTERISTIC ENERGY SCALE AND 
UNIVERSAL SCALING FUNCTIONS 



The exact solution of the Anderson model can be obtained 
using Bethe-Ansatz. Tsvelick and Wiegmann have solved the 
Bethe-Ansatz equations for the Anderson model in an arbi- 
trary magnetic field. Let us summarize their results 14 ^ for 
the symmetric case, which were obtained assuming the con- 
tinuum limit and an infinite bandwidth (th —> oo) for an An- 
derson impurity coupled to a 3D bath of non-interacting elec- 
trons. Since the bath for ISIM is provided by non-interacting 
electrons free to move on a semi-infinite ID tight-binding lat- 
tice, and since we give results for values of t c , td and U which 
are not always small compared to th, one cannot rule out cer- 
tain quantitative differences between the results of Refs. ll4lfl5l 
and our numerical results. This may concern the numerical 
prefactors in the expression of Tk or the constants in the uni- 
versal functions describing the magnetization m(h). How- 
ever, a qualitative agreement should be expected. The Kondo 
temperature of the Anderson model reads 



In the perturbative regime (U < irb%), the conductance is 
well described by HF theory, which yields g = A/(v/t 2 + 
t 2 /v) 2 with v given by Eq. ( f30b . The scale r takes a value 
thf given by Eq. QTT l. 

In the non-perturbative regime (U > nt 2 .), let us revisit 
our numerical results for g using the exact expressions giving 
the magnetization m(h) for an Anderson model (3D bath of 
conduction electrons, wide band limit where th — > oo) which 
is not exactly the Anderson model corresponding to ISIM. 

The conductance g(td) reads 



g{t d ) = sin 2 (2wm(t d )) 



(39) 



where m(td) is the pseudo-magnetization of the nanosystem 
with pseudo-spin 1/2 and reads 



m(t d ) 



(n e (t d )) - (n (t d )} 



(40) 



T, 



K 



(35) 



where F is a prefactor which depends on the definition of Tk 
(which varies 14 from one author to another) and is modified 33 
if the bandwidth of the bath of conduction electrons is finite 
or infinite. F = \/2/tt for the infinite bandwidth Anderson 
model 14 - 15 , while F = Q.ZM^/2jii if the bandwidth is taken 
finite^. 

When U > T (non-perturbative regime), the Bethe-Ansatz 
results for the impurity magnetization m(h) can be divided in 
three characteristic regimes as the magnetic field h increases 
at T = 0. 

A SC regime for low fields (h < Tk) where the magnetiza- 
tion is given by a universal function of h/Tx- 



m(h) = 



2nT K 



(36) 



followed by a LM regime for intermediate fields (Tj. < h < 
\/UT) where m(h) is given by another universal function of 
h I Tk which can be expended as 



m(h) pb - ( 1 l -r— 

2 V ln (£) 



before having a PO regime for strong fields (h > 
noted FO regime in Ref . Il4l) where 



(37) 



m(h) 



1 



2r 



(38) 



becomes independent of the interaction U. 

When U ~ I\ there is direct transition from a non-magnetic 
regime where m(h) « h/T towards the free orbital regime 
where the behavior of m(h) is given by Eq. ( f38l >. 



In the SC limit (fid < Tk) where g and g coincide, this gives 

2 



g(t d ) = sin 2 \td/T k ) ~ f jT" 



9{td) , (41) 



while the conductance g extracted from the NRG spectra reads 

2 / , \ 2 



g(td) 



td/r + r/t d 



(42) 



This yields a relation between the characteristic scale r of 
ISIM and the Kondo temperature Tk of the Anderson model: 



T(U,t 2 c ) = 2T K (U,t 2 c ) 



(43) 



Using for Tk Eq. (f35T > with the finite bandwidth prefactor 
F = 0.364-\/2/7r given in Ref. 33, one gets for t the analyti- 
cal expression 



-(U,t c ) = 0.728 



2t 2 U 



■ exp - 



ttU 
St 2 



(44) 



Eq. d44l describes very well the numerical values of t(U, t c ) 
obtained from the NRG spectra and the condition g(t d = r) = 
1 in the non-perturbative regime, as shown in Fig.[TT]for U > 
t 2 /A. Moreover, the behavior of g calculated for U = 0.25 
and t c = 0.1 is given in Fig. [17] as a function of td when 
t d < t. One can see that the expression g = sin 2 (2td/r) with 
the value of r calculated from Eq. (@4]> (t = 1.582510~ 6 ) 
describes the behavior of g or g calculated using the NRG 
algorithm when t d < r. 

Fig. [18] gives the NRG values of g and g as a function of 
td for U = 0.25 and t c = 0.1 in the LM regime. One can 
see that g s» g around the transmission peak, but becomes 
slightly different when t d . > t. Using Eq. d37l i for m(td), 
the NRG values of g(t d ) shown in Fig.[l8]are not reproduced 
by sin 2 (2irm(td)). We explain this failure by the fact that 
the Anderson model corresponding to Eq. d37l i is not exactly 
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1(T 8 1(T 7 10~ 6 



FIG. 17: (Color online) SC regime: Values of g (red plus), g (green 
cross) obtained from the NRG algorithm and Bethe-Ansatz expres- 
sion sin 2 (2td/r) (solid line) as a function of td for U = 0.25 and 
t c = 0.1. The value of r used in the Bethe-Ansatz expression 
have been obtained from the relation r = 2Tk, with Tk given by 
Eq. HI. 




0.02 



0.01 I — — — — — I 

10~ 7 10~ 6 10~ 5 10~ 4 10" 3 10~ 2 

td 

FIG. 18: (Color online) LM regime: NRG values of g (red plus) and 
g (green cross) as a function of td for U = 0.25 and t c = 0.1. The 
solid line (fit Bethe-Ansatz) corresponds to Eq. ( 145b with the value 
of r = 2Tk given by Eq. J44t and used in Fig. 1171 

the Anderson model corresponding to ISIM. This might give 
different constants in the function given in Eq. d37l >. However, 
we have been able to find a function of X — td/r inspired by 
the form of mQi/Tx) given in Eq. d37l ) and which fits very 
well the values of g in the LM regime: 

^=^(.(2.0175-^^)). (45, 

As shown in Fig. [18] such a fit with the value of t — 2Tk 
used in Fig. [17] allows us to describe g(td)- Eq. ( |45| l gives an 
excellent approximation of the universal curve of g(td/r) in 
the LM regime. 

When td > t c \U (PO regime), g can be described by 
the non-interacting expression A(td/t 2 + t 2 /td)~ 2 , which 
achieves the complete description of g(td, U, tc) in the sym- 
metric case by analytical expressions. 



XV. SUMMARY AND PERSPECTIVE 

When td < r, we have shown that the quantum conduc- 
tance is given by a universal function of the ratio X = td/r of 
two energies. This universal function g(X) = 4(X + X~ 1 )~ 2 
characterizes the non-interacting limit, where the isolated 
nanosystem has two levels of energy Vq ± td with a level 
spacing A = ltd- Those levels have a width V = t 2 when the 
nanosystem is coupled to leads at Ep = 0. Therefore the ratio 
X is also the ratio A/2T. We have found that g(X) remains 
unchanged when the electrons interact inside the nanosystem, 
if one adds a term cx U to the broadening T in the perturbative 
regime. When U becomes larger, there is a non-perturbative 
regime where a more complicated many-body resonance ap- 
pears at Ep. In that case, the relation r = 2Tk which we 
have obtained is consistent with the fact that the width T of 
the nanosystem levels becomes the Kondo temperature Tk- 
We have thus shown that the interaction U leaves the func- 
tion g(X) unchanged in the SC regime, renormalizing only 
the level broadening T. 

To have a conductance which is given by the ratio of 
two characteristic energies and which stays at zero temper- 
ature on a universal curve when this ratio varies is rem- 
iniscent of the scaling theory of localization^!. Recently, 
it has been found 35 using numerical Quantum Monte Carlo 
simulations that the /3(<?)-function characterizes not only the 
non-interacting limit 34 , but also 2D disordered systems with 
Coulomb interactions. We have given another example of a 
universal function which remains unchanged when the elec- 
trons interact, the interaction renormalizing only the charac- 
teristic scale (T for the dimensionless ratio A/T in this study, 
the localization length £ for the dimensionless ratio L/£ in 2D 
disordered systems 3 —). 

When t < td < t c \/U , a local moment is formed in the 
equivalent Anderson model and the universal function g(X) 
is not given by the non-interacting limit. Adapting an analyt- 
ical expression describing the magnetization of the Anderson 
model, we have proposed an analytical form which reproduces 
the universal function g(X) in the LM regime. 

We have shown that the interacting region becomes per- 
fectly transparent when td = t. In the non-perturbative 
regime, this corresponds to a nanosystem level spacing A rj 
Tk- This result shown using ISIM is a particular illustration of 
the minimal realization of the orbital Kondo effect in a quan- 
tum dot with two leads, which has been studied by Silvestrov 
and Imry 36 for more general setups. The role of the conduc- 
tion electron spin is played by the lead index in Ref. 36, while 
it is played by the even and odd orbitals for ISIM. In the two 
cases, the Kondo effect takes place if there are two close levels 
in a dot populated by a single electron, and the conductance at 
T = is zero at the SU(2) symmetric point (td = for ISIM), 
while it reaches the unitary limit G = e 2 /h for some finite 
value Ri Tk of the level splitting A. However, the prediction 
made in Ref. 36, that for temperature T > Tk the conduc- 
tance becomes maximal if the levels are exactly degenerate, 
cannot be valid for ISIM where level degeneracy means no 
coupling between the left and right leads (td = 0). 

We have shown that the quantum conductance of an inter- 
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acting nanosystem coupled to non-interacting ID leads can 
be described with a universal function g(X). This concept, 
with a similar function g(X), must remain valid if one cou- 
ples the nanosystem to non-interacting 2D or 3D leads instead 
of strictly ID leads. This can be understood if one consid- 
ers the Anderson models of pseudo-spin 1/2 particles cor- 
responding to ISIM with leads of dimension D. The Kondo 
physics of such models, where the quantum impurity is cou- 
pled to a bath of dimension D, is qualitatively independent 
of the used bath. A change of the dimension of the leads 
modifies only the dependence of the nanosystem level width 
r upon t c (V oc t 2 c in all dimensions D, but with factors which 
depend on D), and hence the dependence of the Kondo tem- 
perature upon t c . This makes likely that the universal aspects 
of g obtained in a pure ID limit using ISIM do characterize 
also more general spinless models, where the nanosystem and 
the leads would be created in gated 2D semiconductor het- 
erostructures. In that case, ISIM is a simplified model which 
could describe quantum transport of spin polarized electrons 
through an inversion-symmetric double-dot setup, as a func- 
tion of the inter-dot coupling. Since such a coupling can be 



easily varied if the two dots are coupled by a quantum point 
contact, it will be interesting to check whether the quantum 
conductance of such setup is given by a universal function of 
the dimensionless inter-dot coupling td/Tx when T — > 0. For 
observing the orbital Kondo regime using such a setup, a large 
capacitive inter-dot coupling will be necessary. It will be also 
interesting to introduce the spin 1/2 of electrons in ISIM for 
studying the role of f<j upon the SU(4)-Kondo effect, as we 
have studied its role upon the SU(2)-Kondo effect using spin- 
less fermions. The possibility of observing SU(4)-symmetric 
Fermi liquid state in a symmetric double quantum dot system 
with strong capacitive inter-dot coupling has been discussed 
inRef.H3 

Eventually, this study was restricted to the symmetric case, 
leaving to a following work the study of the asymmetric case 
(Ep ^ 0, Vg 7^ —U/2), where the role of td upon the 
valence-fluctuation fixed point remains to be investigated. 

We thank Denis Ullmo for very useful discussions and the 
RTRA "Triangle de la Physique" of Palaiseau-Orsay-Saclay 
for financial support. 
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